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LTN/Inlets and Nozzles Branch 
Overview 


Computational and Experimental 
Research for Advanced Inlet and 
Nozzle Concepts 


*-Design & Analysis Tools 
*Fundamental Flow Physics Investigations 
Small & Large Scale Component Testing 
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Turbulent CFD Validation Experiments 


High-order large eddy simulation of shear 
layers, mixing layers and separated flows _ v.nasa.gov 
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Primary Experimental Facilities 
Inlets and Nozzles Branch 


¢ CE-22: Nozzle Static Test Facility 

¢ (CW-13, CW-17): Free Jet Facilities 

¢ 1°x1’ Supersonic Wind Tunnel (SWT) 

¢ 15cm x 15cm (and new 17-cm 
axisymmetric) SWT 

¢ AeroAcoustic Propulsion Laboratory 


¢ 10’x10’ Supersonic Wind Tunnel 

¢ 9’x15’ Low Speed Wind Tunnel 
” ¢ 8’ x6’ Supersonic Wind Tunnel 

Sipe ¢ Propulsion Systems Laboratory 
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NAS Supercomputing 
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BLI Propulsor Technology Genesis and Applications 


x _ Credit: United Technologies Research Center / NASA 
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Boundary Layer Ingesting Inlet and Distortion Tolerant Fan (BLI2DTF) 


Robust Design Experiment Layout S 
Boundary Layer Ingesting Inlet 
NASA 8’x6’ Transonic Wind Tunnel 
ot 
oo s » 


Convergent / Divergent Bump Ce BP 
To Provide BLI Shape Factor B 


OBJECTIVES: 

¢ Design an experiment which accurately simulates the integrated 
propulsion system Boundary Layer Ingesting Inlet and Distortion 
Tolerant Fan stage on a candidate Hybrid Wing Body/ Blended 
Wing Body vehicle. 


Universal Propulsion 
Simulator Drive Rig 


¢ Validate that the performance and operability of the UTRC 
provided boundary layer ingesting inlet and distortion tolerant fan 
stage meets the design requirements & specified performance 
metrics. 


Raised Test 
Section Floor 


Boundary Layer 
ae 7 F . : ee = Fast-Acting, Flow Calibrated 
* Utilize the experimental database obtained during this test activity Variable Area Nozzle 
for CFD validation and code development (i.e. both aerodynamics 


and aeromechanics/structural dynamics codes). 


Distortion-Tolerant Fan Stage 
+ Rotating AIP Rake Array 

+ Distortion-Tolerant Fan Stage 

+ Rotating Fan Exit Rake Array 


¢ Establish the aeromechanic behavior of the Distortion Tolerant 
Fan (DTF) and compare to pre-test predictions. 
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ACCOMPLISHMENTS: Raised Floor 


, - Diffuser V: 
* Very successful wind tunnel test -operability, research goals ee ear 


met 

¢ Unique new BLI test bed infrastructure successfully 
demonstrated 

¢ High quality data acquired —virtually no instrumentation 
failures 

¢ Mach number range: 0.125 and 0.500 to 0.783 

¢ Boundary layer thickness range: 3.65 to 5.00 inches 

¢ Fan stage pressure ratio range: Windmill to 1.37 

e Research results to be reported at the 2017 AIAA Propulsion 
and Energy Forum, Atlanta, GA and 2018 AIAA SciTech 


Pre-Entry Diffusion Ramp 


Credit: NASA/ David Arend 


POINT OF CONTACT: 
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The TBCC Propulsion Concept sal 


¢ The Turbine-Based Combined Cycle (TBCC ) concept involves an over/under 
configuration with a split in the propulsion flowpath. 
— From take-off to Mach 4, a turbine engine provides propulsion. 
— Above Mach 4, a dual-mode supersonic combustion ramjet /scramjet (DMSJ) 
engine provides propulsion. 
¢ Inlet Mode Transition is the procedure by which the flowpath split is achieved 
and controlled as the turbine engine shuts down and the dual-mode 
ramjet/scramjet engine takes over. 


Mach < 4 


Mach > 4 
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Turbine-Based Combined Cycle Engine Inlet and Controls vasa 
research in the GRC 10’x10’ SWT 


RESEARCH GOALS 


Proof of concept of over/under split flow inlet 


- Develop performance & operability 
database for both the turbine and scramjet 
inlets 


- Demonstrate mode transition and control 


Low-speed ramp 
& flowpath 


‘ A  High- speed ramp 


& flowpath 
Mach 


¢ Validate CFD predictions for both flowpaths nurbine Engine 
¢ Develop realistic distortion characteristics 
throughout the mode transition Mach number 
range Dual-Mode 


« Testbed for mode transition controls research 
& for integrated inlet/engine propulsion 


systems TEST APPROACH - 4 Phases 


Inlet performance and operability characterization 


Completed . System Identification of inlet dynamics for controls 
. Demonstrate Control strategies for smooth & 


stable mode transition without inlet unstart 


Add engines/ nozzle for integrated system test 
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Axisymmetric Shock-Wave/Boundary Layer 
Interaction Experiment 


Objectives: 


. Obtain mean and turbulence quantities through a M=2.5 
SWBLI of sufficient quantity and quality to be considered 
as a CFD validation dataset. Initial efforts will focus on a 
Mach 2.5 2-D (in the mean) interaction with follow-on 
efforts investigating 3-D interactions. Both attached and 
separated interactions will be considered. 

Technical Challenge: 

° Very few Shock Wave/Boundary-Layer Interaction 

(SWBLI) experiments reported in the open literature 

meet the rigorous criteria required to be considered as a 

CFD validation dataset. This is particularly true for 

experiments with detailed turbulence measurements. 


Deliverables/Products: 

"Complete characterization of an axisymmetric 
shock-wave/boundary-layer interaction 

" High quality CFD validation dataset including 
experimental uncertainty quantification 

= New axisymmetric supersonic tunnel facility, 
including PIV measurement system 


Approach/Key Milestones: 

¥ Completed flow surveys using conventional pressure 
probes 

¥ Completed testing using dynamic surface shear 
stress film (S3F) and fast response Pressure 
Sensitive Paint (PSP) in collaboration with Innovative 
Scientific Solutions, Incorporated (ISSI) 

¥ Completed flow surveys using normal hot-wire 
anemometry 

o Complete flow surveys using slant hot-wire 
anemometry 

o Complete flow surveys using particle image 
velocimetry (PIV) 


POINT OF CONTACT: 
LTN/Dave Davis, david.o.davis@nasa.gov 


17cm Axi-Supersonic Wind Tunnel (in ERB W6B cell) 
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OBJECTIVES 


Perform validation experiments to collect turbulent flow 
measurements in the region of an injection hole and at several 
locations downstream. 


Implement advanced non-intrusive diagnostics technique for 
measuring temperature fluctuations in a plane with high spatial 
resolution and high temporal response rates 


Data will be used to validate CFD methods (RANS, LES, DNS), to 
develop boundary conditions for film cooling, to develop more 
accurate models for turbulent heat flux, and to improve 


TECHNICAL CHALLENGES 


Accurate measurement of fluctuating 3-component velocities and 
temperature at high speed and high temperature close to a surface 
and/or over an injection hole 

Temporal and spatial resolution capabilities of measurement 
techniques 


APPROACHI/KEY MILESTONES: 3-Steps 
v STEP 1: Perform a flat plate test examining single hole and multiple 


hole patterns using SW-6 facility. Low speed (up to Mach 0.2) room- 
temperature crossflow with heated film injection flow (AT up to 100 
degs F). 


STEP 2: Collect fluctuating temperature data from a heated 
convergent jet (up to 1300 degs F) tested on the Small Hot Jet 
Acoustic Rig using a new advanced RAMAN spectroscopy capability 


STEP 3: Perform a flat plate test installed in a hot free-jet (either 
downstream of a nozzle in the SHJAR or the NATR) with ambient 
temperature film cooling flow injected through a single or multiple 
holes using a new advanced RAMAN spectroscopy capability. 


Complementary CFD simulations for each step at pre-test to help 
define the experiment and post-test to validate/improve codes and 
models 


-Transformational Tools & Technologies (T*) Project 


Turbulent Heat Flux (THX) Validation Experiment 


PRODUCTS 


¢ High-quality flow field data including: Mean velocities and 
temperatures, kinetic energy, Reynolds stresses, surface 


temperature, dissipation rate, turbulent viscosity or length scales, and 


turbulent heat flux 


« New advanced RAMAN spectroscopy to provide improved non- 
intrusive diagnostic capability at GRC for measuring mean and 
fluctuating temperatures in a plane at high spatial resolution 


PARTNERS/CUSTOMERS: Boeing, P&W, AFRL 


Step 3-Nozzle Geometry 


0.5” Diameter Injector, 30° inclined 
to axial direction. Center of hole 
Is 1.5 in from nozzle exit. 


A 


Plate is 12” by 2.68” by 0.1”. 
Thickness will likely need to 
increase. Capability for supporting 
thermocouples. 


POINTS OF CONTACT 
CFD; Nick Georgiadis, 3-3958, georgiadis@nasa.gov 
RAMAN/PIV: Mark Wernet, 3-3752, mark.p.wernet@nasa.gov 
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INCREASING FIDELITY 


Tools and Methods vasa 


Inlets and Nozzles Branch 


Method of Characteristics (MOC) design codes: Several in-house computer programs for 
designing high performance supersonic inlets and nozzles. 


Inlet and Nozzle Performance Analysis Code (IPAC & NPAC) - predicts performance for a given 
geometric design, can be used to design preliminary inlet & nozzle systems and to make 
subsequent performance analyses. 


SUPersonic Inlet design and analysis tool (SUPIN)- uses analytical, numerical, and empirical 
methods to design the geometry and estimate the aerodynamic performance of supersonic inlets. 
Steps include: 
— establish the parametric design of the inlet 
— efficiently model the geometry and 
— generate the grid for CFD analysis with design changes to those parameters. 
RANS Codes: 
¢ Wind-US - a general purpose, applied CFD tool. It solves the Navier-Stokes equations of 
fluid mechanics for steady state and unsteady flows, along with supporting equation 
governing turbulent and chemically reacting flows. NASA GRC personnel partner with the 
USAF AEDC and Boeing to develop, support, and validate the code. 
¢ FUNS3D — a general purpose, applied CFD. It solves the Navier-Stokes equations of fluid 
mechanics for steady state and unsteady flows, along with supporting equations governing 
turbulent flows. Developed and maintained by NASA Langley. 
LES and Hybrid RANS/LES (Wave-Resolving LES (WRLES) & Glenn Flux Reconstruction 
(GFR codes)) — high-order codes for propulsion flow modeling. 


www.nasa.gov 


What is Wind-US and NPARC? sal 


The NPARC (National Project for Applications-Oriented Research in CFD) 

Alliance is a partnership of NASA GRC and USAF AEDC, with significant 

participation of the Boeing Company, to provide the state-of-the-art in proven 

CFD technology. 

¢ The NPARC Alliance has been in existence for 20 years; significant evolution 
from PARC code in early 1990s to Wind-US Version 4.0. 

¢ The primary product of the NPARC Alliance is the Wind-US (US=unstructured) 
code with accompanying utilities and extensive on-line documentation. It is 
made available to all U.S. organizations approved to receive export-controlled 
information. 

¢ Wind-US is a general-purpose, multi-zone, compressible Navier-Stokes CFD 

code. It is designed to provide the state-of-the-art in proven CFD code 

technology. 


It is the workhorse flow solver for NASA GRC Inlets and Nozzles branch. 


¢ Wind-US code distinguishing features: very detailed documentation and 
validation suite, broad range of physical models. 
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Validation of Combustion, 


Structured & Unstructured Solvers 
Mach 4.5 Boundary Layer: 


u/U 


1.2 


0.8 


0.6 


0.4 


0.2 


BL profile at: 
x=1.79 ft, Re,=3.81 million 


Str-Str, SST 
Str-Str, SST w/Comp Corr, Press Dilat OFF 


Str-Str, S-A Burrows—Kurkov — Stagnation Temperature (K) 


Str-Uns, SST 

Str-Uns, SST w/Comp Corr, Press Dilat OFF 
Str-Uns, S-A 

Uns-Uns, SST 

Uns-Uns, S-A 

Coles 23: Fence Trip 

Coles 22: Fence Trip 

Coles 61: Sand Trip 

Coles 62: Air Jet 

Experiment 
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Dutton-Goebel Supersonic Mixing Layer: 


Struct. 


Un.-St. 


Burrows-Kurkov Reacting Mixing Layer: 


y (mm) 


L ° Goebel Data 
10 SST, structured 
[| sseseeeeee SST, unstructured 
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Wind-US Validation Website 


+ NASA Portal 


Validation Web Site Search: 


ae 
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[Home] Archive [Resources] [History] [Contacts] 


NPARC Alliance Verification and Validation Archive 


Last Updated: Thursday, 22-Jan-2015 13:02:29 EST 


e 
The archive consists of cases. Each case represents a geometric configuration and flow condition (i.e. ONERA M6 wing, MADIC 3D ( } t t d ft ft bs 
nozzle). A case can be classified as either a verification or validation case depending on whether the comparison data for the case is ns YUuUC ure je excerp e 
exact (analytical) or experimental, respectively. Each case contains one or more studies. Each study consists of one or more runs of 
the Wind-US code (or another CFD code). Throughout some of the cases and studies are detailed examples that provide a 


step-by-step description or tutorial on using WIND-US and associated tools. Table 3. Unstructured, 3-D grids for unstructured solver 
; Name Azimuthal Width (deg) | Azimuthal Cells | Centerline Scheme | Total Cells 
Downloading Files from the Archive 
Str-Uns_90deg-v1 90 20 triangular-prisms 1,006,847 
The cases of the archive include: Str-Uns_90deg-v2 90 20 tiangular-prisms 1,957,622 


Mach Grid Gas Description 


0.00 SG IG Sod’s Shock Tube : 

0.05 SG IG Incompressible Driven Cavity Tags: 

0.10 SG IG Basis incompressinte Laminer Flat Plate SG = structured grid 
0.13 SG IG Driver-Seegmiller Incompressible Backward-Facing Step OG= orid 
0.15 SG IG Fraser Subsonic Conical Diffuser UG = unstructured grid 
0.20 SG IG Incompressible, Buice Axisymmetric Diffuser IG = ideal gas 

0.20 OG IG NLR Airfoil with Flap CH= 


0.20 SG,UG IG Incompressible, Turbulent Flat Plate 


0.20 SG IG Laminar Flow over a Circular Cylinder 

0.20 SG IG Ejector Nozzle 

0.21 SG IG Low-Subsonic S-Duct 

0.30 SG IG Steady, Inviscid Flow in a Converging-Diverging Verfication (CDV) Nozzle 

0.30 SG IG Subsonic Annular Duct 7 : Ne De 
040 SG IG Square Jet Injection Va) 
046 SG,UG IG Sajben Transonic Converging-Diverging Diffuser 

0.73 SG IG RAE 2822 Transonic Airfoil 

0.80 SG IG MADIC 2D Axisymmetric CD Boattail Nozzle 

0.80 SG IG MADIC 3D CD Boattail Nozzle 

0.84 SG IG Transonic, ONERA M6 Wing 

0.97 SG,UG IG Acoustic Reference Nozzle with Mach 0.97, Unheated Jet Flow 


1.30 SG IG Normal shock at Mach 1.3 
1.82 SG,UG CH _ Hydrogen-Air Combustion in a Channel 
2.00 SG,UG IG Mach 2.0 Flow over a 15-Degree Wedge 


a.) Full computational grid. Zoom view of nozzle region. i *"c.) Zoom view of nozzle exit. 


2.00 SG,UG IG Seiner Nozzle with Mach 2.0, Heated Jet Flow Lid 
222 sG IG Supersonic Axisymmetric "submerged" Jet Flow Figure 4. 3-D, unstructured grid for Acoustic Reference Nozzle, grid Uns-Uns_90deg-v2. 
2.35 SG IG Conical shock on a 10 degree cone at Mach 2.35 


2.44 SG CH Burrows and Kurkov Supersonic Mixing/Combustion 
2.50 SG,UG IG Oblique shock on a 15 degree wedge at Mach 2.5 

250 SG,UG IG Prandtl-Meyer 15 Degree Expansion Corner at Mach 2.5 
450 SG,UG IG Mach 4.5 Flow over a Flat Plate 

UG IG Mach 5.0 Shock Boundary Layer Interaction 

7.00 SG IG Mach 7, Laminar 15-degree Ramp 

15.00 SG IG Hypersonic Cylinder 
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Propulsion Airframe Integration CFD 


Potential need for additional validation data and CFD tool refinements 
Cruise Pt, 


¢« NASA analysis of Quiet Supersonic Mach 1.42, AOA 1.7 
Technology (QUeSST) design was performed: 
— Large differences (>5%) in NASA FUN3D 
vs. LM CFD++ inlet recovery and 
distortion. 
— Differences attributed to predicted flow- 
field on top of vehicle. 


¢ At GRC, the QueSST project marks the first 


time full vehicle aerodynamic analysis was 92% Recovery 97% Recovery 
performed using FUN3D on a mature vehicle Fun3D CFD++ 
design. 


¢ The QueSST aircraft is unique in that it 
contains an aft-fuselage, top-mounted inlet 
installation, with complex shock-wave 
boundary-layer interactions that are sensitive to 
geometry and are difficult to predict. 


¢ Depending on the success of this first test, 
more investment/validation data sets likely 
needed to bring PAI best practices up to the 
same maturity level as sonic boom analysis. 
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HWB Propulsion Airframe Integration CFD 


* Developed a design approach for hybrid “ 2 -—— sx 
wing body (HWB) integrated airframe of ov .. “=e, ‘ch 
and propulsion by parameterizing the Mi v \ eee | 
planform and airfoils at control sections |; “Ta SN 
of the wingbody ES tae =) — : = 
* Used Go-Flow code, a 3D unstructured "SS . ——— 
N-S flow solver, performed Euler ee - — Sy 
analysis’; flow through propulsion 3 : - ——_—_— 
represented by a simple block model oer | _ 
‘ 0.814286 Baseline Optimized 
and a proper static back pressure at the 0 842867 | 
fan face Dams: 
- Developed adjoint-based optimizer to 044295 
improve aerodynamic performance 02 | _ 
¢ Reduction in induced and wave drags 2 ) | | 
were achieved, resulting in ~10 counts 4.2 oe 
reduction in total drag under the oS te ait + vs as = mat 1 | \ 
trimmed condition | | ~ “pot rN 
* Demonstrates mutual interference : = ne, 
between airframe and propulsion a a = 
system is significant and can be ——— . la ee OE a a 
accounted for using this computational . = 
procedure 


Optimized geometry comparison: Baseline (left, blue) vs 
Optimized (right, red) of the integrated configuration 
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LES Research for Propulsion 


Objectives: 


«Develop tools and techniques for simulating propulsion flows 
using large-eddy simulation (LES) and hybrid Reynolds- 


Deliverables/Products: 
« High quality LES databases for model development 


Averaged Navier-Stokes/LES methods. » 2D mixing layers 
"Create detailed databases of shear layer and shock- " Hot jets 
wave/boundary-layer interaction flows to guide new model = Shock-wave/boundary-layer interactions 


development « Diffusers, separated flows 


= Accurate turbulent inflow boundary conditions and 


Technical Challenge: best practices for simulations 


While LES is a conceptually simple technique there are 
numerous issues that must be resolved before it becomes 
routine analysis tool; inflow turbulence specification, wall 
modeling, hybrid RANS/LES models, numerical scheme, etc. 


Approach/Key Milestones: 
¥ Incorporate inflow turbulence generation methods 


v¥ Develop detailed flow field solutions for mixing 
layers 


¥ Develop detailed flow field solutions for hot jets 


o Develop detailed flow field solutions of shock- 
wave/boundary-layer interactions 


Incorporate wall models 


Develop detailed flow field solutions of diffusers, 
separated flows 


© ©) 
0 5 10 x 16 20 250 5 10 x 15 20 25. 


Density contours of mixing layer (close-up view of splitter tip) 
for high-order LES (left) and low-order RANS (right) 
[from AIAA2015-2939, Mankbadi & Georgiadis] 


POINT OF CONTACT: 
LTN/Jim DeBonis 


-Transformational Tools & Technologies (T?) Project www.nasa.gov 


National Aeronautics and Space Administration 


Informing Turbulence Models Using ILES 


PROBLEM 

Accurate prediction of shock wave/boundary-layer interaction (SBLI) 
remains a challenge using the current industry-standard turbulence 
models. This is largely due to the lack of understanding of SBLI 
characteristics like low-frequency motion of the shock-separation 
system, three-dimensionality, and corner separation where the 
Boussinesq approximation starts becoming suspect. 


OBJECTIVES 

Using the insights provided by Implicit LES (ILES) in the physics of 
low-frequency shock-separation oscillations, three-dimensionality, 
and corner separation, inform the current turbulence models to gain 
improved predictions and/or show the need for development of new 
higher order turbulence models. 


APPROACH 

Highly-resolved meshes will be used in SBLI simulations performed 
with ILES to, first, acquire statistical quantities to study the unsteady 
characteristic of an SBLI and second, to calculate budget of exact 
equations of turbulence kinetic energy (TKE), rate of dissipation, 
specific rate of dissipation, and stress transport and compare with 
their RANS counterparts. An 8-degree SBLI exhibiting evidence of 
three-dimensional shock-separation system will be investigated. 
Effect of corner separation will be explored by including a sidewall in 
the simulation. The code, FDL3DI, employs bandwidth and order 
optimized WENO with Roe scheme for inviscid flux and a sixth-order 
compact scheme for the viscous fluxes. The time integration will be 
performed using implicit Beam-Warming scheme. Digital filtering 
method will be used to generate inflow turbulent boundary layer and 
and the shock is imposed using the Rankine-Hugoniot relations by 
specifying pre- and post-shock conditions. 


POC: Manan Vyas (GRC) 


DETAILS OF PRELIMINARY WORK 

Preliminary simulations using a flatplate have been completed and verified the 
suitability of turbulent boundary-layer generation for the present flow conditions 
and calculation of components of the budget terms like production, transport, 
dissipation, dilatation, etc. However, the SWBLI length is smaller in comparison 
to the experiment and future work is aimed at uncovering the source of the 
discrepancy. 


SIGNIFICANCE 

This two-pronged approach will provide physical insights on the role that 
unsteadiness and three-dimensionality play in SBLI and lead to improvements in 
our current models and spur the development of new models that incorporate 
the relevant flow physics. 


Streamwise Velocity 
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¢ Separation is present in the ILES, but in the near-wall region 
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Glenn Flux Reconstruction (GFR) Navier-Stokes Code 


PROBLEM 


Current industrial CFD codes used for solving aeropropulsion flows with 
complex geometries generally use unstructured finite-volume methods and are 
typically no higher than second-order accurate. These methods constrain the 


accuracy and efficiency of current CFD codes. 


OBJECTIVES 
At Glenn, develop a high-order 3D CFD code based on the Flux 


Reconstruction (GFR) method for solving the Navier-Stokes equations on 
brick element meshes to demonstrate accuracy and efficiency improvements 


over standard methods. 
APPROACH 


The GFR code provides a simple, efficient, and easy to implement method 
for solving fluid flow problems on mixed element unstructured grids and is 


accurate to an arbitrary order. The GFR code is to be used for solving 
aeropropulsion flows with complex geometries. This code has been 


expanded from 2D inviscid flows to 3D viscous flows on brick elements. 


While extending the 2D code to 3D, additional infrastructure has been 


embedded within the GFR code to easily allow for future additions such as 


3D hybrid meshes and non-uniform P-adaptive meshes. 


RESULTS 

The isentropic Euler vortex problem was computed to verify the 
accuracy of the FR method for inviscid flows. Results of the 
study found the expected P+1 order accuracy through one 
period and super-accuracy through 100 periods demonstrating 
that vortical flow structures can be accurately predicted by GFR 
for long amounts of time with minimal unwanted numerical 
dissipation introduced into the system. 


SIGNIFICANCE 

This 3D Navier-Stokes GFR code uses the FR method to 
efficiently compute high-order solutions on brick element 
meshes. Its capabilities are being extended to allow for hybrid 
element meshes and more generalized aeropropulsion-oriented 
flows. Results indicate that the GFR code appears to be a viable 
candidate for accurate ILES of aerodynamic flows for complex 
geometries. 


POC: Seth Spiegel (GRC) 
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Aggregate Wall Time / DoF 


Comparison of the aggregate wall times per DoF versus 
time step for each of the vortex simulations. This shows 
that for each time step, the computational cost per DoF is 
essentially independent of P: 


www.nasa.gov 


